Interacting black holes on the brane: the seeding of binaries 
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We consider the evolution of sub-horizon-sized black holes which are formed during the high energy 
phase of the braneworld scenario. These black holes are long-lived due to modified evaporation and 
accretion of radiation during the radiation dominated era. We argue that an initial mass difference 
between any two neighbouring black holes is always amplified because of their exchange of energy 
with the surrounding radiation. We present a scheme of binary formation based on mass differences 
suggesting that such a scenario could lead to binaries with observable signatures. 
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Cosmology in the braneworld scenario has inspired 
much activity in recent times. The feasibility of a large 
extra spatial dimension is the cornerstone of the Randall- 
Sundrum (RSII) braneworld modeljlj in which all the 
standard model fields are confined to our observable 3- 
brane, except for gravity which can propagate also in the 
bulk. The braneworld description of our universe entails 
an early high energy phase during which the evolution of 
the universe is significantly altered. A particular feature 
of interest in RSII cosmology is the evolution of primor- 
dial 5-dimensional black holes. Black holes forming out of 
the collapse of horizon-sized density perturbations dur- 
ing the high energy phase obey a different evaporation 
law0 compared to 4— dimensional black holes. It has 
been shown recently that accretion from the surrounding 
radiation bath can dominate over evaporation during the 
high energy phase This leads to prolonged survival 
of these primordial black holes with multifarious cosmo- 
logical and astrophysical consequences. There exists the 
possibility that these black holes could be a significant 
fraction of cold dark matter. Several observational con- 
straints on primordial black holes get modified in the 
braneworld scenario, as has been shown recently Q. 

Black holes as dark compact objects in galaxy haloes 
are the target of many recent and ongoing observa- 
tions||- A large number of these black holes are expected 
to be in the form of binaries. Several interesting obser- 
vational evidences ranging from lensing effects to hyper- 
velocity star ejections have been put forward to bolster 
claims in support of black hole binaries 6] in a wide mass 
spectrum. Recently, it has been proposed that sub-lunar 
mass binaries comprised of braneworld black holes could 
produce detectable gravitational waves in their coalesc- 
ing stage0- Gravitational lensing experiments do indeed 
leave open the possibility of the existence of sub-lunar 
mass black holes in certain mass ranges 0. A forma- 
tion mechanism for binaries in standard cosmology based 
on the inhomogencous spatial distribution of primordial 



black holes has been proposed[2j- 

In this Letter we present a scheme of binary forma- 
tion for primordial braneworld black holes based on mass 
differences. We consider sub-horizon sized black holes 
in the RSII model formed during the high energy phase 
whose evolution has been studied in Refs. 0, El • We shall 
consider the interaction of two neighbouring black holes 
which arc initially well separated for their gravitational 
attraction to be approximated by Newtonian dynamics. 
The physical distance between two such neighbours in- 
creases with the Hubble expansion of the high energy 
radiation dominated era. These black holes exchange en- 
ergy via the processes of evaporation and accretion with 
the radiation bath, and through it, with each other. The 
initial mass ratio of two such black holes is on average 
proportional to the square of the ratio of their forma- 
tion timesQ. We will see from the evolution that the 
initial mass difference between two such black holes al- 
ways increases during the radiation dominated era. We 
will then argue that such mass differences will facilitate 
the formation of binaries during the standard low energy 
phase via three-body gravitational interactions. We will 
consider an example of sub-lunar mass binary0 forma- 
tion through this scheme. 

Let us first briefly state some essential results of 
Refs. HQ. 

The mass of a primordial braneworld black 
hole formed at time to with initial mass mo during the 
high energy phase (t < t c with t c = 1/2) grows due to ac- 
cretion (evaporation is negligible, except for black holes 
with very small initial mass, or very low accretion effi- 
ciency /, i.e., except for (Af /^4) 2 < 2A/(2B-l)(t /t c ), 
where A ~ 3/((16) 3 tt), and B = 2//tt, with < / < 1) 
such that 
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This equation holds up to t ~ t c assuming radiation 
domination persists, i.e., the radiation density pr — 
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(3m|)/[327r(t c + t)t] stays greater than the sum of the 
energy density in all black holes, pbh- The formation 
time is related to the initial mass by 



(2) 



with the 4-dimensional Planck mass (time) labelled as 
m 4 (£4). Assuming accretion to be effective only during 
the high energy phase, these 5-dimensional black holes 
evaporate as 
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during the standard low energy regime t > i c , with m max 
being the maximum mass a black hole reaches via accre- 
tion. [Recently, there has been a conjecture of rapid evap- 
oration through conformal modes for large braneworld 
black holes with size greater then the AdS radius Z [lo| . 
However, our focus will be on much smaller black holes.] 
Now let us consider the evolution of two neighbouring 
black holes, BH1 and BH2, having masses mi and 7712 
respectively, with initial masses mi < m,2 , or t\ a < ti a . 
BH2 forms at a physical distance do from BH1. Their 
physical separation d grows at a rate oc t 1 / 4 . We shall 
assume that do > I, and that the gravitational potential 
4> = [m/(m|d)][l + {2l 2 )/(3d 2 )} < 1 so that exchange 
of energy between them via gravitational waves can be 
neglected during the high energy regime. The change 
of mass of each black hole is effected by the net sum of 
its Hawking evaporation depending on its temperature, 
and the accretion of radiation depending on the mean 
radiation density in the universe 3j. Moreover, due to 
the presence of BH2, the net rate of change of the mass 
?7i 1 of BH1 gets a further contribution proportional to the 
product of 777,2 and the solid angle subtended by BH2 on 
BH1, and vice- versa. In other words, BH1 feels the local 
difference from the mean radiation density as a result 
of an accreting (or evaporating) BH2, and similarly for 
BH2. Thus the evolution equation for each of them can 
be written as 
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for i, j = 1,2 [(n/h) = (8/37r) 1 /2(//;4 ) i/2 (m . /m4) i/2 
is the Schwarzschild radius for 5-dimensional black 
holesQ]. The first and second terms on the r.h.s rep- 
resent accretion from the average radiation density, and 
evaporation respectively The last term arises from 
the local inhomogeneity in radiation density due to the 
j-th black hole. Considering both BH1 and BH2 to be in 
their accreting phases at the formation time to 2 of BH2, 
it is possible to see from Eqs.(@J that for BH1 the first 
term is comparable to the interaction (third) term when 
{d/r 2 ) ~ {[t-h }/[t-t 2o }) 1/2 . Thus, for a while after its 



formation, BH2 is able to suppress the growth of BH1. If 
on the other hand, BH2 forms at a time when BH1 begins 
to evaporate, with or without the aid of the interaction, 
BH2 registers enhanced growth due to the locally denser 
radiation coming from the evaporating BH1. 

The effect of two interacting black holes on their evolu- 
tion can be more precisely formulated as follows. Defin- 
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where 



\t/t c , Eqs.(@J can be written as 
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can be regarded as a coupling constant between BH1 and 
BH2 dependent on their initial separation do- Now, in 
terms of the variables s = ln(t/to) and Xi = m.i/t 1 / 2 , 
Eq.Q can be cast in an autonomous form, i.e., 
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where accents denote derivatives with respect to s. In the 
following we confine our attention to small enough values 
of <7, such that g 2 x\x 2 < 1. This is indeed the condition 
for Eq.JTJ to be well-behaved. For two equal mass black 
holes (xi — x 2 = x), this condition becomes g < g+ — 
l/x. This limit coupling strength in fact corresponds to 
the limit of coalescence, since the black hole radius r is of 
the order of the separation d for g = g+. Below the critical 
value g c = [2{2B — l) 3 /27] 1 / 2 of the coupling constant, 
the dynamics of two equal mass black holes has a stable 
fixed point x f ■ Our interest is in two unequal mass black 
holes which are initially well separated for gravitational 
interaction to be subdominant. If one considers a small 
initial mass difference, then the two black hole masses 
will evolve as x\ = x[s) — e(s) and x 2 — x(s) + e(s), with 
e(s) e(0)exp (J Q A(x(u))duj , where the instantaneous 
Lyapunov exponent 
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is larger than the constant (2B — l)/2. It follows that 
any small initial mass difference blows up exponentially 
in the first stages of dynamical evolution. In particular, 
the symmetric fixed point x f is always linearly unstable 
against a perturbation of the form 8x\ = —Sx 2 . 

It is thus clear that the physical solutions of Eqs.JSJ) 
exhibit disequilibration for black holes masses, i.e., if for 
two neighbouring black holes, mi (to) < m2(to), then 
mi(t) <§; m2 (t) for t > to. Requiring the era t with a 
number density hbh of black holes of average mass m(t) 
to be radiation dominated, i.e., pbh = Jtiubh < PR, the 
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mean strength of the coupling g between two such black 
holes should be such that 
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In general however, there could be spatial inhomo- 
geneities in the formation of black holes. The growth 
rate of any given black hole is hence impacted by the 
mass spectrum and spatial distribution of all neighbour- 
ing black holes. For a system of two black holes BH1 and 
BH2 formed at eras t\ and t 2o with the initial mass dif- 
ference (A 12 )o = Tn,2 — mx , and when both are in their 
accreting phases, the mass difference at time t 3> i 2o is 
given by (from Eqs.(l) and (2)) 



m 2o 




m 2o 



(A u )(i) - (A 12 ) > 

1 - B "\-( 1 . 



(10) 



with the equality sign holding for g = 0, i.e., when the 
interaction term can be neglected. 

For the present analysis, we are assuming that black 
hole accretion is restricted to the high energy phase. 
Hence, when the universe exits the high energy phase, 
the masses of BH1 and BH2 are frozen to mi and m 2 
(say). Hawking evaporation continues to be negligible 
for a long time beyond t c [jj]. Consider a spherical region 
of radius d in which the average matter density (con- 
tributed by BH1 and BH2) is given by 
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Since radiation density in this region falls off as a~ 4 , the 
sphere of radius d will become matter dominated at a 
time tf given from Eqs.(l) and (2) by 
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For t > tf the matter dominated region with radius d gets 
cut off from the background radiation dominated expan- 
sion which continues up to the era of matter-radiation 
equality t eq 3> t c . For t > tf, the black holes BH1 
and BH2 will form a bound systemQ. Assuming that 
a bound system (or binary) gets formed only in the low 
energy regime (tf > t c ), it follows from Ea. (|12|l that 
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With the initial density fraction a mo (= Pbh{U)) / PB,{to)) 
for black holes formed with mass mo, the mean initial 
separation is do = (37r/8) 1 / 2 (ro/am 3 ). [Two black holes 
with initial separation less than that given in Eq. I|13|) will 
form a bound system in the high energy phase. Accretion 



of radiation will be negligible after binary formation, with 
a resultant smaller lifetime. In the present analysis we 
shall not consider details of such a scenario.] The size of 
this bound system df is related to the initial black hole 
separation by df — do(i c Ao) 1 ^ 4 (^/Ac) 1 ^ 2 - Using Eqs.(l), 
(2) and 113, 
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with its mean value given by df = (A/oiml)(to/t c ) 1+B l. 

Binary formation proceeds as a result of 3-body in- 
teraction. Let us choose three neighbouring black holes 
to have mass diffrences at t c given by to 2 — mi = Ai 2 
and m,3 — m 2 = A 2 3. Such a scenario will arise generi- 
cally since the initial mass differences are due to horizon 
sized collapse at different eras, and possible interactions 
through the radiation background can only amplify such 
mass differences between neighbouring black holes during 
the high energy regime. BH2 will become gravitationally 
bound to BH3 at time tf with an elliptical orbit whose 
major axis 



a f 



(15) 



The role of BH1 is to provide the tidal force required to 
prevent head-on collision between BH2 and BH3. The 
minor axis 6/ of the binary is proportional to the tidal 
force times the squared free fall time0, i.e., 
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where dp is the perpendicular distance of BH1 from 
the mid-point of the BH2-BH3 axis. Although BH2 
binds to BH3 at time tf, mass differences hamper the 
possibility of BH1 being gravitationally bound to BH3. 
Even if the three black hole separations are equal be- 
fore binary formation, the total region encompassing the 
binary and BH1 will become matter dominated at a 
time (t md )i 3 much later than tf, given by (i TO d)i3 > 
*/[( TO 3 + TO 2„)/(mi +m 2o )] 11/i ~ B , (since using Eqs.(2) 
and (11), tf oc n / 4 ) with the equality sign signifying 
g = 0. But during the interval tf < t < (i m( z)i3, BH1 
continues to move away from the BH2-BH3 binary with 
the Hubble flow, i.e., dp registers growth by a factor 
[{t m d) 13 /tf] 1 / 2 . The gravitational interaction between 
the binary and BH1 at t = (t m( i)i3 is much weakened 
compared to the strength of the BH2-BH3 interaction. 

The eccentricity of the binary orbit can be obtained 
from Eqs. i|15|) and l|16fl. For uniformly distributed black 
holes formed around the era t with masses ~ m , the 
mean eccentricity is given by 



1 - r; 2 2 



m 2 - A12 
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1/2 



(17) 
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where -q is an O(l) geometrical factor|9j. The Newtonian 
approximation that we have used for describing the gravi- 
tational interaction between black holes forming a binary 
is valid if the gravitational potential \<p\ <C 1, which for 
the BH2-BH3 binary reads (mi/mf) <C df(l — e). Using 
Eqs.(l), (|14|l and (|17|l one can verify that the Newtonian 

analysis holds good for (^/io) 25-1 *C rfo^J 3 [{m,2 — 
A 12 )/2m 2 + A 23 )] 2 . 

As a particular example of the scenario outlined by 
us, consider three black holes formed about the era 
to ~ 10 26 t4 with average mass mo ~ 10 24 m4, radius 
r - 10 27 / 4 and A - O(mi). We have chosen i = 10 3 %, 
about the maximum allowed by present experimental 
bounds. Setting the accretion factor B = 1/2, one gets 
m max ~ m(t c ) = 10 26 TO4 ~ 10 21 g. The closure bound for 
a population of black holes with mo ~ 10 24 m 4 can be ob- 
tained using the prescription of Clancy et al. |4j , and turns 
out to be a mo < 6 x 1CT 18 for the initial energy fraction. 
Choosing the maximum allowed value of a mo the mean 
initial separation is do ~ 10 6 ro ~ 1cm. Nevertheless, 
the validity of the Newtonian approximation \(f>\ <C 1 is 
readily ensured for a choice do ~ 10 3 ro for the specific 
system of three black holes we are considering for binary 
formation. For these parameters, binary formation takes 
place at tf ~ 10 12 t c , with the major axis aj ~ 10 7 cm. 
Note here that we have estimated the mean values for 
the above parameters assuming a uniform distribution of 
black holes. Inhomogcncitics in the spatial distribution 
for black holes would change these values on a case-to- 
case basis. However, at least some of the binaries would 
have parameters close to the mean values. These would 
then fall in the category of sub-lunar compact objects dis- 
cussed in Ref. , and gravitational waves emitted during 
their coalescence are within the observational range of 
third generation gravitational wave detectors. 

To summarize, we have studied the effect of interaction 
of horizon-sized black holes formed during the radiation 
dominated high energy phase of the RSII braneworld sce- 
nario. We have seen that the processes of black hole 
evaporation and accretion mediated via the radiation 
bath causes the mass difference between two neighbour- 
ing black holes to be always amplified during the high en- 
ergy era. We have presented a scheme by which binary 
formation takes place through three-body gravitational 
interaction during the standard low energy era. This 
scheme of binary formation based on mass differences 
could be viewed as an additional possible mechanism to 
the scheme based on spatial inhomogeneitiesQ- In a re- 
alistic scenario, both spatial and mass inhomogencitics 
are expected in any given configuration of black holes, 
and these separate effects will interfere to either abet or 
inhibit binary formation depending upon the particular 
configuration. The specific example we have furnished 
suggests that the parameters for such binaries might lie 



within the range of sub-lunar compact objects amenable 
for detection by forthcoming observations of gravitational 
waves such as by the EURO detector Q. It is important 
to note that we have considered black hole accretion only 
during the high energy phase. If accretion of radiation 
continues during the standard low energy regime up to 
t cq , as claimed in Refs.ll], it might be worthwhile to 
investigate the possible formation of more massive bina- 
ries through primordial black holes. One should keep in 
mind that black holes also grow by accreting mass from 
the surrounding matter distribution in the matter domi- 
nated era. Finally, the possibility of accreting the homo- 
geneous and presently dominating cosmological or scalar 
field energy 12], if realized, could push up the masses fur- 
ther. 
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